Bounds in proton-proton elastic scattering at low 



o 
o 
o 

(N 
•*-> 
o 

O 

(N 



> 
^1- 



O 



momentum transfer 

A. T. Bates and N. H. Buttimore 

School of Mathematics, University of Dublin, 

Trinity College, Dublin 2, Ireland 
atbates@maths . ted . ie , nhb@maths . ted . ie 

February 1, 2008 



Abstract 

o, 

We present a bound on the imaginary part of the single helicity-flip amplitude for 
spin 1/2-spin 1/2 scattering at small momentum transfer. The variational method 
of Lagrange multipliers is employed to optimize the single-flip amplitude using the 
values of otot, °"ci and diffraction slope as equality constraints in addition to the 
inequality constraints resulting from unitarity. Such bounds provide important in- 
formation related to the determination of polarization of a proton beam. In the 



^ . case of elastic proton collisions the analyzing power at small scattering angles offers 

a method of measuring the polarization of a proton beam, the accuracy of the po- 
larization measurement depending on the single helicity-flip amplitude. The bound 
obtained on the imaginary part of the single helicity-flip amplitude indicates that 
the analyzing power for proton-proton collisions in the Coulomb nuclear interference 
region should take positive nonzero values at high energies. 



1 Introduction 

The proton spin puzzle has intrigued experimentalists and theorists since the surprising 

result from the EMC experiment at CERN in 1988, which found a smaller than expected 
*A.T.B. and N.H.B. are grateful to Enterprise Ireland for partial support under Scientific Research 
Project SC/96/778 and International Programmes IC/1999/075, IC/2000/027. 
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contribution to the spin of the proton from the component quarks. The question, "where 
does the spin of the proton come from?" remains unanswered [|l], || . Recent data suggests 
a value of ~ 20 — 30% for the fraction of the spin carried by the up, down and strange 
quarks. The contribution from the gluons and from the orbital angular momentum of 
the quarks and gluons is not completely known. The Relativistic Heavy Ion Collider 
at Brookhaven National Laboratory plans to probe the proton structure using the deep 
inelastic scattering of protons at high center-of-mass energies (y/s = 50 — 500 GeV) and 
momentum transfers (pr > 10 GeV/c) @. To measure the contribution of the gluons to 
the spin of the proton, with sufficient accuracy, a polarized proton beam with a beam 
polarization error of 5% is necessary. 

One method of measuring the polarization of a proton beam uses the analyzing power in 
elastic proton collisions at small scattering angles. The analyzing power An for a proton 
and the transverse single spin asymmetry A are related through the expression 

A N P = A (1.1) 

where P is the beam polarization; for 100% beam polarization the asymmetry and ana- 
lyzing power are equal. The beam polarization can be measured by counting the scatters 
with the beam polarized up (N^) and then down (N^) in a polarimeter with a known 
analyzing power An'- 
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(1.2) 

A N 



W + iV-l- 

The analyzing power An expressed in terms of the s-channel helicity amplitudes is 

djG 4:71 

A 's = -Kr^) lm[4tl * i+ *' + * , ~ M - (L3) 

The reduced ratio, r 5 , of the helicity single-flip to imaginary non-flip amplitude is defined 

as 

r 5 = ^=x-^-. 1.4 
' t Im 0+ 



The analyzing power An for the CNI region can be written as follows, when the 
transverse total cross section spin difference is neglected [HI: 



N -t (k - 2Imr 5 )(t c /t) + 2(plmr 5 - Rer 5 ) 
N ~ m (t c /ty - 2(p + 5){t c /t) + 1 + p 2 + /? 2 1 ' j 

where t c = — (87r«)/cr tot , 5 is the Coulomb phase, k + 1 = ji = 2.7928 is the magnetic 
moment of the proton, and (3 2 represents the forward hadronic double spinflip contribu- 
tions expected to be negligible at high energies. Apart from the photon pole term, the 



t-dependences of helicity nonflip and flip electromagnetic and hadronic amplitudes due 
to form factor and nuclear slope effects are not expected to play a significant role in the 
amplitude ratios featuring in the asymmetry. Inside the CNI region interference between 
the electromagnetic and hadronic amplitudes is most prominent. An important contribu- 
tion to the maximum of An, in the CNI region (|t| < |t c |), comes from Imr5 in the form 
of n — 1 — 2Imr 5 . At larger momentum transfers outside the CNI region (|t| > |t c |) the 
analyzing power, containing hadronic terms only, is essentially 



A V^t 2 Qlmr 5 - Rer 5 ) 

An = 1, 2 ; m — L6 

m 1 + p 2 + p z 

and because the value of p is expected to be no more than about 10% at RHIC energies, 
the analyzing power provides a good indication of the value of Rer 5 just outside the 
CNI region. The contribution of Re rs to the maximum of the analyzing power would be 
reasonably well known. A bound on Imrs which satisfies fi — 1 — 2 Imrs > ensures that 
the maximum analyzing power in the CNI region is positive. 

The optimization technique of Lagrange multipliers, extended by Einhorn and Blanken- 
becler [|, [| [7| to include equality and inequality constraints in the context of scattering 
theory, is used to derive the bound 

1/2 

Otot 



,, Cm^l^-l) x/,f/-) (1.7) 



where 

(l+ 9 t(2 + 9gt/8)) 1/2 

h[t) ~ (TT^) ' (L8) 

The partial wave expansions for the observables, a tot , a e \, and the hadronic slope pa- 
rameter g to be defined in Section 2, are included as equality constraints in the system, 
unitarity is presented as partial wave inequality constraints and the modified imaginary 
single-flip amplitude, the function to be optimized, is input as the objective function. 
Before proceeding to the derivation of the bound in Section 3 we first introduce the s- 
channel helicity amplitudes, expanded as partial wave series in the helicity representation. 
Errors on the derived bounds are discussed in Section 4 followed by concluding remarks 
in Section 5. 
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2 Helicity Amplitudes and Observables 

For the elastic scattering of two protons at CM energy ^/s and CM momentum k 



s — 4m 2 /2, there are sixteen helicity amplitudes in general, each a function of s and t. 
The number of independent amplitudes is reduced under the following relations || 

Parity conservation 

(a;a 2 | |AiA 2 ) = (-if~ A (-a; - a' 2 | 1 - A a - a 2 ) (2.1) 

Time reversal invariance 

(AiA^I |A!A 2 ) = (-lf~ A (AiA 2 | |A;A 2 ) (2.2) 

Identical particle scattering 

(A'iA 2 | lAiAa) = (-1) A ^ (AgAil |AaAi) (2.3) 

where A = Ai — A 2 , \i = \[ — A' 2 . 

The sixteen helicity amplitudes reduce to two non helicity-flip amplitudes 0i and 3 , 
two double helicity-flip amplitudes 2 and 04, and one single helicity-flip amplitude 05, 
with partial wave expansions || || : 

X (5, t) = (+ + |0| + +) = ^E (2J + 1) (Zo J (*) + /£(*)) 4) (6) (2.4) 
03 (a, t) = (+ - 101 + -) = ^ E (2J + 1) (//(*) + f 2 J 2 (s)) d J u (6) (2.5) 



2 (s, t) = <+ + |0| - -) = ^ E (2 J + 1) (/£( a ) - tf( a )) 4 (0) (2.6) 
04 (3, *) = (+- |0| -+> = t^E (2J + 1) (/£(«) - (0) (2.7) 

5 ( S , = (+ + |0| + -) = ^ E (2^ + 1) /£(«) d{ (6) (2.8) 



where // (i = 0, 1, 11, 22, 21) denote s-channel partial wave amplitudes, Imff = a/, 
Re// = 6/ and 2 = cos# = 1 + t/2k 2 . In the Coulomb Nuclear Interference (CNI) 
region, t « — 0.0012(GeV/c) 2 , it is convenient to express the five helicity amplitudes in 
terms of Jacobi polynomials. To define the d J x (8) function in terms of Jacobi polynomials 
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it is suitable to separate the space of A and \i into four regions A, B, C, D [I(J. In region 
A, where A + /i > and A — /i > 0, the relation is 



(J + A)!(J-A)! fl + z\^~ (l-z\—*~ 



xPtx^\z), (2-9) 



\\ (J + ii)\(J-ii)\ V 2 J \ 2 

and J — A = 0, 1,2, ■ • -. Equivalent forms in the other regions are obtained by use of 

symmetry relations |], . 

In region B, where A + /i > and A — ji < 0, use 



<W = (-1) A ^<aW- 
In region C, where A + < and A — < 0, use 

d{,{9) = (-1) A -" <£ A _ M (0) . 
In region D, where A + /i < and A — /i > 0, use 



(2.10) 



(2.11) 



(2.12) 



Expressing the d J X(l {6) functions in terms of Jacobi polynomials the five independent he- 



licity amplitudes can be written as 

M*,*) = ||E(2^+l)(/o J (^) + /n(^)) ^i° ,0) W 
= ^E(2^+l)(/i 7 i(^)-/o J (^))^ '° ) W 

4 (s,t) = 



Ak 


L 




A, 




v^Vl 


-z 2 



(0,2), 



,0), 



4A; 



£(2.7+1) 



7 "ViW^W 



J 



(2.13) 
(2.14) 
(2.15) 
(2.16) 

(2.17) 



where z = cos# = 1 + t/(2k ). The Jacobi polynomials have the properties |TT 



T(a + n + 1) 
T(a + l)n! 



^ P^Hz) = o-m r ( m + n + a + / j + 1 ) D (a+m,/3+m) / \ 



dz 



T(n + a + (3 + l) 



(2.18) 
(2.19) 
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2.1 Observables 

In proton-proton elastic scattering the spin observables can be written in terms of the 
five helicity amplitudes, <f>i, . . . , $5. The observables are vital to the optimization, since 
each observable can be included as an equality constraint in the optimized system. In 
the derivation of the bound we use three equality constraints, a to t, cr e \ and g, and two 
inequality constraints related to unitarity. 

2.1.1 Total Cross Section 

The first equality constraint uses the total cross section. The optical theorem, 

Im<f> + (s,t)\ t=0 = ^a tot (s), (2.20) 

is used to express the total cross section as a partial wave expansion given by 

^(s) = ^(2J + l){(a J +a J u ) P^°\l) + (a{ + a J 22 ) (2.21) 

where Im (j)+(s, t) = (Im <fii(s, t) +Im 3 (s, t))/2 is the imaginary spin average helicity non- 
flip amplitude. Using property ( ETTB ) the normalized dimensionless total cross section can 



be expressed partial wave expansion: 



^o = ^E( 2J + 1) K (s) + a{ (s) + o&( a ) + a J 22 (s)} (2.22) 
where A Q = (k 2 /n) a tot . 

2.1.2 Slope of the Imaginary Non-Flip Amplitude 

The slope of the imaginary non-flip amplitude has been used in bounds for other spin 
dependent elastic collisions [0. We find it convenient to use the imaginary part of the 
spin-averaged amplitude at a particular value of t. The second equality constraint employs 
the imaginary spin average non-flip amplitude at a particular t inside the Coulomb Nuclear 
Interference region, written as a Taylor expansion: 



Im0 + (s, t) w Im 0) +t ( -^Trn0 + (s, t) 



(2.23) 

t=o 



where \t\ is sufficiently small so that inclusion of the linear term in the Taylor expansion 
is an accurate approximation. Use of properties (|2.18| ) and (|2.19| ) leads to the partial 
wave expansion for the imaginary non-flip amplitude: 

Im Ms, t) = ^J2(2J+ 1) {a J +a{ + a J u + a J 22 } (l - | J(J + 1)^ (2.24) 



where ( = —t/k 2 . The logarithmic derivative of the imaginary spin average non-flip 
amplitude, 



9 = Jt Xn Im< M*'*)l«=° = Im^Q) {j t lmMs ^ 



t=o 



can be expressed as 



-^Im0 + (s,t) 



g Im0 + (s, 0) . 



(2.25) 



(2.26) 



t=o 



The Taylor expansion for Im0 + (s,t), given by Equation ( |2.23|) , can thus be written as 



]m<f> + (s,t) =Im0+(s,O) {!-(-*<?)} 



(2.27) 



2.1.3 Elastic Cross Section 

The third equality relates to the elastic cross section, expressed as a partial wave expansion 
by integrating the differential cross section over momentum transfer t: 



0"ei(sj 



o dt da<M) 
-4k 2 dt 



(2.28) 



Expressing the differential cross section in terms of helicity amplitudes allows us to write 
the elastic cross section as 



7T 



Using the expression 



2k 2 s 



f° rft{|0i| 2 + |0 2 | 2 + |03| 2 + |04| 2 + 4|0 5 | 2 } • (2.29) 

J-Ak 2 1 } 



t = -2k 2 (l - z) 



(2.30) 



the t variable can be replaced with the z variable, where t is the momentum transfer, k 
is the center-of-mass three-momentum and z = cos 6. The elastic cross section expressed 
as an integral over z becomes 



>ii 2 +i0 2 i 2 +i03i 2 +i</>4r 



+ 4|0 5 | 2 } 



(2-31) 



To express the elastic cross section as a partial wave expansion, the integrals 

' +1 dz\<j>i{s,t)\ 2 



are calculated, where i — 1, • • • , 5. The integration formula [13 



(1 - z) a {\ + zf P^\z) Pt P \z) dz 



m 7^ n 

2 a+/3+l r(q+n+l) T(/3+n+l) r 

n\ (a+/3+2n+l) V(a+f3+n+l) ° m ' n ' 



(2.32) 



may be used to find the integral 

r+l 



dUe)<$J0)dz 



>M, J 



(2.33) 



2J + 1 

leading to a partial wave expansion for the normalized dimensionless elastic cross section, 
defined as E e i = (k 2 /n) cr cl : 



s el («) = y: (2j + 1) {|/ f + 1 A J f + |/iif + |/ 2 J 2 | 2 + 2 |/ 2 i| 2 } 



(2.34) 



2.2 Imaginary Single-Flip Amplitude 

The imaginary single helicity-flip amplitude, modified by a kinematical factor, is the 
objective function in the system. Before optimization we must first express the single-flip 
amplitude in a suitable form. The imaginary amplitude, 



(2.35) 



Im0 5 (M) = |^£(2J+ 1)^4,(0) 



written in terms of Jacobi polynomials and normalized by a kinematical factor, becomes 

Im0 5 (s,t) y/s. 



Im i 



(1 - z 2 )V2 4k 



£(2J + 1) 



(2.36) 



In the CNI region, the Jacobi polynomial P^i(z) expanded as a Taylor series is 



(2.37) 



2 = 1 



where z — 1 + t/(2k 2 ) and £ = —t/k 2 . Using properties ( 2.18 ) and ( |2.19|) , the Taylor 
series for Pj 1 ^^) about z — 1 is 



F£?(*)« J(1-MJ(J+1)-2]1 , 



c 



(2.38) 



and thus, inside the CNI region 



I m 5 ^£(2J+l) ^±Ij^l-| [J(J +l)-2]j a 2 J x 



At small collision angles, 



??? 



3 ^(l-z 2 ) 1 / 2 ' 
the ratio Imrs = mIm05/(^/— t Im0 + ) is 



(2.39) 



(2.40) 



Im r§ 



m Im( 



k lm(f> + (s,t) 
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2.3 Unitarity 

The partial wave amplitudes obey the following unitarity inequalities [14 



u( 


= a J 


- 1/oT > 




ui 


= 4 


- \fi\ 2 > o 




vl 


= 4x 


— \fii 1 2 — fii 


2 > 




= a 22 


~ I/22I — I/2I 


2 > 



(2.41) 
(2.42) 
(2.43) 
(2.44) 



where f£ 



0,1,11,22,21) denote the s-channel partial wave amplitudes, Im// 



Cl; 



and Re// = bj. Combining Equation ( p.41j ) with Equation (|2.42|) leads to the inequality 
X J = U{ + Ui where 

X J = al + al-\fi\*-\fi\ 2 >U (2.45) 



and the inequality W J = V/ + V/ follows from the combination of Equation (p.43|) and 



Equation ( 2.44 ) where 



W 



«i 1 +^2-iAir-i/ 2 j 2 r-2i/ 2 j 1 r>o. 



J |2 



cJ |2 



;2.46) 



For the elastic scattering of spin on spin 1/2 particles there are two independent 
helicity amplitudes, a flip and a non-flip amplitude, with partial wave expansions whose 
partial wave amplitudes obey unitarity relations similar to relations ( 2.41Q and ( 2.42 ). The 
unitarity relations (|2.43|) and (|2.44j) are characteristic of spin 1/2 - spin 1/2 scattering, 
the fi x term coming from the single helicity- flip amplitude ^5. 



3 Optimization 

Equipped with partial wave expansions for the observables and partial wave inequality 
relations, representing unitarity, we are in a position to optimize the modified helicity 
single-flip amplitude Im0 5 . We follow the variational technique of Einhorn and Blanken- 
becler § by constructing a Lagrangian consisting of an objective function and a set of 
equality and inequality constraints. We use the full set of constraints, cr to t, crci, 9 and 
unitarity, although a bound on Im0 5 , with fewer constraints, can be derived fl~5|1 . 
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3.1 Lagrange Formalism 



The normalized dimensionless total cross section A , expressed as an equality constraint, 
is included in the Lagrange function along with the normalized dimensionless elastic cross 
section E e i, written as an equality constraint, the imaginary spin average non helicity-flip 
amplitude Im0+(s, t) at a fixed small \t\ value, also expressed as an equality constraint, 
and the partial wave unitarity relations written as inequality constraints. The modified 
single helicity-flip amplitude Im 05 is introduced as the objective function in the Lagrange 
function: 



C 



Im 5 + a 

+/3 



eJ |2 
'21 1 



+ 7 



Aq-Y, (2 J + 1) {a J +a{ + a J n + a J 22 ) 
j 

s el - £(2j + 1) (l AT + IAT + 1 Ail 2 + I Ail 2 + 2 1 a 

J 

Im + - ^ £ (2 J + 1) {a J Q +a{ + a J n + a J 22 } (l - ^ J(J + 1) 

+ £(2J + K + < ~ I Ail 2 - I Ail 2 - 2|/ 2 i 12 
j 

+ Y,(2J+l)\j K + aM/oT-lAT) 
j 

where a, /3 and 7 are equality multipliers. The inequality multipliers, Aj and /xj, are by 
definition non-negative and £ = —t/(k 2 ). In the high energy or large J limit only the 
leading order J terms are included and the Lagrange function of Equation fl3.1|) becomes 



(3.1) 



C 



Im 05 + a 

+/3 
+7 



A - 2 ^ J [4 + af + a J n + a J 22 ) 
j 

E el - 2 £ J (l/of + |//| 2 + l/^l 2 + |// 2 | 2 + 2|/ : 
./ 

Im 0+- £ J{a^ + af + 0^ + 0^} ^1 - | J 2 ^) 

j 

+2^JA J (a J + af-|/ J | 2 -|/ { 



J |2 
21 1 



fj |2 
'22 I 



2|Ail 2 



and 



Iuti 



(3.2) 



(3.3) 



The system is optimized by taking first and second derivatives with respect to the real 
and imaginary partial wave amplitudes, 6/ and a/. This gives the optimized set of partial 
waves, at some fixed t in the CNI region; 

hi 



Mi. 



(3.4) 
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01 l + 2Aj V ; 

°«=^= — W — (3 - 6) 

and 

« 2 J i = ^ f 1 " I ^ (3.7) 

where Aj = Aj/2/3, /ij = /ij/2(3, r\ = —a/ (2(3), r 2 = —7/(4/3) and /3 > for a maximum 
(or (3 < for a minimum). 

3.2 Unitarity Classes 

Optimization under the four constraints imposes the following condition: 

6f = 0Vi=>// = of + 6f = of, (3.8) 

that is, there is no contribution from the real partial wave amplitudes. The imaginary 
partial wave amplitudes therefore obey the following unitarity inequalities: 

X J = a J Q - a J G 2 > (3.9) 

W J = a J n - a J u 2 - a J 21 2 > (3.10) 

where and ad = a J 22 - It is natural to divide the partial waves into two classes, 

one with contributions from the interior unitarity class I and the other with contributions 
from the boundary unitarity class B. For the X J unitarity inequality the interior and 
boundary unitarity classes are defined as 

I x = {j\Xj>0, Aj = 0} , B x = {j\Xj = 0, Aj>0} . (3.11) 

Likewise for the W J unitarity inequality the interior and boundary unitarity classes are 

I w = {J\Wj>0, jlj = 0} , B w = {J\Wj = 0, jlj>0} . (3.12) 

3.2.1 I x and B x Unitarity Classes 

The interior unitarity class, 

I x = {j\X J = a J -ai 2 >0, Aj = 0} , (3.13) 
11 



under the four constraints, is expressed as 

I x = {j|0 < a J Q < 1, Aj = 0} . (3.14) 

Equation ( |3.5|) with A j set to zero enables us to write the imaginary partial wave amplitude 

j 
• 



Oq, in the interior unitarity class, as 



<4 = r 1 + r 2 \l-^J 2 j . (3.15) 

The constraint < Oq < 1 restricts the values of the equality multipliers, T\ and r 2 , to 
< T\ + r 2 < 1 and r 2 > 0. The number of partial waves J are thus limited to 



0<^<-^l + -j. (3.16) 

The boundary unitarity class B x splits into two sub-classes, B x ° and : 

— ► 5 X ° = {j|a^ = 0, Aj > 0} 

fi x = { J I Xj = a J - a J 2 = 0, Aj > o} 



= {j\a^ = l, Aj>0} 



(3.17) 



In the boundary unitarity class B x ° the imaginary partial wave amplitude is equal to 
zero and from Equation (|3.5| ) the inequality multiplier Aj is given by 

Aj = -(r 1 +r 2 )+r 2 ^ J 2 >0. (3.18) 

The 5^° class begins at J 2 = Mf = 4/( (l+n/r 2 ), and for J > Afi+1, with < ri+r 2 < 
1 and r 2 > 0, the inequality multiplier Aj is positive. Therefore the boundary unitarity 
class B x ° is non-empty for J > Mi + 1 but with Oq = 0, for all J in this unitarity class, 
there are no contributions to the observables from this unitarity class. The imaginary 
partial wave amplitude Oq is equal to unity in the boundary unitarity class B Xl and from 



Equation (|3.5|) the inequality multiplier Aj is given by 



Aj = (ri + r 2 )-l-r 2 ^J 2 . (3.19) 
By definition Aj > and the value of J, in the boundary unitarity class, is limited to 

J2 £ 4 pL+^i) ,3.20) 

but with < r\ + r 2 < 1 and r 2 > 0, J 2 is negative, or J is complex and therefore the 
boundary unitarity class B Xl is empty. 
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In summary, the unitarity classes, I x and B x °, are non-empty and the unitarity class 
B Xl is empty; 

I x = {j|0 < a J < 1, < J< Mi}, (3.21) 
B x ° = [j\a J =0, Mi + 1 < J < M 2 ) (3.22) 



where Mi = Floor ^JA/C (1 + r 1 /r 2 ) , M 2 = Floor ^8/( and ( = -t/k 2 . The 
Floor [x] function gives the greatest integer less than or equal to x. 

3.2.2 I w and B w Unitarity Classes 

The interior unitarity class I w under the optimization becomes 

I w = { J | W J = a J u - a J u 2 - a J 21 2 > , jlj = o} , (3.23) 

Substituting Equations ( |3.6|) and (|3.7|) , with jlj = 0, into the interior constraint a{ x — 
a ii 2 ~ a ii 2 > leads to the equation; 

f 2 (J)=d 1 + h 2 J 2 + a 3 J 4 + a i J 6 >0 (3.24) 

where a x = (r 1 +r 2 ) (l_( ri +r 2 )), 5 2 = r 2 ( (2(r 1 + r 2 )-l)/4-l/(64/3 2 ), fi 3 = C/(256/3 2 )- 
r|£ 2 /16, a4 = — ( 2 /(64(3) 2 , and only positive J solutions are allowed. The solution is of 
the form 

7 2 „ 2 4 / 1 , r A 2 ,,^2 



< J 2 < ^ - [1 + - ^ M 2 X (3.25) 
where r/ 2 ~ 1. The boundary unitarity class B w is written as 

B w = [j\W J = a J n - a J n 2 - a{ 2 = 0, fa > o} . (3.26) 

The constraint af x — a{ 2 — a J 2 2 = can be written as a quadratic equation: 

+ jij + f 2 (J) = (3.27) 



where 



The solutions are 



fz(J) = a 1 + a 2 J 2 + 5 3 J 4 + 5 4 J 6 . (3.28) 



1 



A^-|±^l-4/ 2 (J)-l| . (3.29) 

The function / 2 (J) is negative for J > M x = Floor[4/£ (l + ri/r 2 )] and consequently 
fij is positive for such J values. By definition jlj > 0, therefore the positive solution is 
chosen; 

»J = \ {V 1 - 4 /2(J)-l} • (3.30) 
13 



To summarize, both the classes, I and B w , are non-empty: 



J 2 



7 J 2 



B 



w 



{J\t 
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J 2 
l ll 



; J 2 



> 0, < J < Mi} , 
0, Mi + 1 < J < M 2 } , 



(3.31) 
(3.32) 



with r? 2 = 1, where Mi = Floor (1 + n/r 2 ) , M 2 = Floor and C = -t/A; 2 . 

It is important to notice that with 7] 2 = 1 both interior unitarity classes, I w and I x , are 
non-empty over the same region, J e [0, Mi]. Similarly the boundary unitarity classes, 
and S , are non-empty over the same region, Mi + 1 < J < M 2 . In other words 
there is no mixing of unitarity classes, all classes are either interior unitarity classes, 
I = I x U I w , or boundary unitarity classes, B = B x U B w , for a given J. 



3.3 Solution of Interior Unitarity Class 



Consider the set of interior unitarity classes, I = I U I . The inequality multipliers, Xj 
and jlj, in the interior region are equal to zero. The imaginary partial wave amplitudes 
are therefore written as 



ri + r 2 1 



and 



4 



^ ^ C T 2 N 



1 



(3.33) 

(3.34) 
is the maximum 



= 0, 1, 11, 22, with < J < Mi, where Mi = Floor J4/( (1 + n/r 2 ) 
J in the interior unitarity class. In this case the contributions to the observables and 
to the objective function Im0 5 solely come from the interior unitarity class /; Aq = A , 

= Im0 + , = S el and Im^g 
section is reconstructed by substituting Equation 

Ah 



lm(j)l_ = Im (/>_i_, Egj = S c i and Im0g = Im0 5 . The normalized dimensionless total cross 

into 



A = ^2 J ( a o + a ( + a n + 



'22 



to give 



j=o 



Mi 

A = 8J2-J 

j=o 



r l + r 2 1 



C j2 



(3.35) 



(3.36) 
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The Euler-Maclaurin expansion [I6| for large J is used to write the normalized dimen- 



sionless total cross section Aq as an integration over J: 



A Q « 8 ( Ml dJ\ (n + r 2 ) J-rS J 3 ) (3.37) 



M 2 



4 

{8{n + r 2 )-r 2 CMl} . (3.38) 



2 

Similarly the imaginary spin average helicity non-flip amplitude Im <f>£ is reconstructed 
by substituting Equation ( |3.33j ) into 

Mi 



Im0 + = £ J {4 + a{ + a J n + a J 22 } ( f - C -J 2 ) (3.39) 



j=o 



4 



to give 



r Mi ( t C 2 1 

Im0+ « 4/ dJ{(r 1 +r 2 )J-(2r 2 + r 1 )^J 3 + r 2 }-J 5 \ 
Jo [ 4 lb J 



(3.40) 



« M 2 |2(n + r 2 ) - (2r 2 + r^M 2 + ^^i 4 | ■ (3-41) 

The dimensionless normalized elastic cross section E^, by substituting Equations ( |3.33|) 
and ( g3g ) into 



Mi 

E el = 2 ]T J (a£ 2 + of 2 + a J u 2 + a J 22 2 + 2a J 21 2 ) , (3.42) 



j=o 



is reconstructed: 



E el « |4(r 1 + r 2 ) 2 M 2 -(n + r 2 )r 2 CM 1 4 + ^M 1 6 | 

M 4 C <f (T 2 1 

+ — ^7 U - - M, 2 + — Mi 4 } (3.43) 
64/3 2 \ 6 1 128 1 J V ; 

The modified imaginary single-flip amplitude Im0 5 is reconstructed by substituting Equa- 
tion fiCHD into 



leading to 



Im0 5 = £./ 2 f 1 "!^) 4i (3-44) 

W 5 ^f/(l-^f. (3-45) 
For large J the modified imaginary single-flip amplitude is written as 
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An expression for the equality multiplier (3 is found by solving Equation ( |3.43| ): 

P = 12 — — 1 m ■ ( 3 - 47 ) 

8 {S cl - (4(n + r 2 ) 2 M2 - (n + r 2 )r 2 C M x 4 + Mf) } 

Rewriting the modified imaginary single-flip amplitude one obtains, 

Im0 5 « ^-^(ri + ra^-Cri + r^raCM^+^M^j 7 

X ^L (i _ £ M a + il M 4 1/2 ( 3 48) 

4 I 6 1 128 1 I V ; 



The equality multipliers, r± and r 2 , are found by solving Equations ( p.38| ) and (|3.41|) . The 
solutions are given by 

= AU (l-SIm^M.) 
36(1 -2Im0+/A o ) 



and 



r 2 = ^ ? (3.50) 

72(l-2Im0 + /A)) 2 



where £ = —t/k . The equality multiplier /3, with solutions for n and r 2 , is expressed as 



^ _ 9(A - 2 Im cj> + ) yjl - 2 Im + /A o + 36 Im ^/A§ 
2A Cx/72E cl - 2Ag C/(l - 2 Im0 + /A o ) 



The optimized modified imaginary single-flip amplitude, expressed as a function of r±, r 2 
and (3, becomes 



_ ~ (A o -2Im0 + ) v /l/2-2(l-Im0 + /A o ) Im</> + 

Im0 5 = — — Y = (3.52) 

4 ^o C J36E cl — A§C/(1 — 2 Im + /A o ) 



with 



12/ 9 Im0 + \ 
-'max = — 11-2 — J . (3.53) 



For low momentum transfers the imaginary spin average non-flip amplitude Im0 + , ex- 
panded to order t, is written as 

Im0 + «^(l+#t) . (3.54) 

Under this approximation the maximum J inside the interior unitarity class is independent 
of t and in the limit t — > 0, the number of partial waves is finite where 



Jm &x = Jl2gk. (3.55) 
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The equality multipliers in the low t limit become 



n - ^ (±+M) , ,3.56) 



72g 2 k 2 V t 

r2 = ~r$m (3 ' 57) 



and 



/72S cl - 2A 2 /(gk 2 ) 

The upper bound on |Imr 5 |, where |Imr 5 | = m | Im0 5 |/(/c Im0 + ), can be expressed 



analytically: 



\[2m kg 

Imr5 < 



A 2 

18E el 5- x (3.59) 



A V^ cl 2^ 2 



where 

(l+^(2 + 9^/8)) 1/2 
*(*) = . (3.60) 

The variable is finite at £ = and changes 'slowly' over the CNI region. Writing 
A = k 2 a tot /i{ and E el = k 2 <T e \/ir, enables the bound on | Imr 5 | to be expressed as 

1/2 

at 



\ m ,, | <mj-g\ 367r f ffcl - 1 ] x . (3.61) 



'tot 



3.4 Results 

The value of the bound on | Im r§ | is given in Tables [l] and @ with the values of the 
equality multipliers. The most noticeable feature of the bound is its size at low momentum 
transfers, having a value of 0.89 at yfs = 52.8 GeV, t = —0.001 (GeV/c) 2 . The optimized 
partial waves, at s/s = 52.8 GeV and t = -0.001 (GeV/c) 2 , is shown in Figure 0. The 
partial wave series terminates at J = 231 which is the upper J limit, M\, for the interior 
unitarity class /. When considering both the interior and boundary unitarity classes, 
values of J > Mi are permitted. 

The bound on | Imr 5 |, under the approximation 

'"afe- (3 ' 62) 

with momentum transfers in the CNI region is expressed as 

|Imr 5 | < m J| x h(t) (3.63) 
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Table 1 : Im r 5 


max jngidg the interior 


region at t 


= -0.001 (GeV/c) 2 


Vs (GeV) 


r i 


To 


P 




1 ImrJ 

| ' O | 


19.4 


-12.54 


12.77 


90 


81 


0.97 


23.5 


-12.56 


12.79 


117 


98 


0.92 


30.7 


-12.02 


12.24 


158 


131 


0.92 


44.7 


-11.22 


11.44 


217 


195 


1.05 


52.8 


-11.53 


11.76 


293 


231 


0.89 


62.5 


-11.56 


11.79 


358 


276 


0.86 



Table 2: | lmr 5 \ max inside the interior region at t — —0.01 (GeV/c) 2 . 



(GeV) 


n 


T2 


P 


^max 


Im rs 


19.4 


-1.25 


1.49 


85 


81 


0.97 


23.5 


-1.05 


1.27 


111 


98 


0.91 


30.7 


-1.00 


1.22 


150 


131 


0.92 


44.7 


-0.92 


1.14 


204 


195 


1.05 


52.8 


-0.95 


1.17 


276 


231 


0.89 


62.5 


-0.94 


1.18 


337 


276 


0.86 



and in the zero momentum transfer limit, t — > 0, the bound on | Imrs| is finite and can 
be expressed analytically as 

I Im rs | < m . (3.64) 

This approximation generates a 'stricter' bound on |Imrs|. The results are given in 
Table El. 



0.1 r 




50 100 150 200 50 100 150 200 



Figure 1: a{ (k = 0,1,11,22) and a J 2X in the interior class; yfs = 52.8 GeV, t = 
-0.001 (GeV/c) 2 . 
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Table 3: | lmr 5 \ max , with an approximation for g, over the CNI region. 



^ (GeV) t = (GeV/c) 2 t = -0.001 (GeV/c) 2 t = -0.01 (GeV/c) 2 



19.4 


0.803 


0.805 


0.825 


23.5 


0.805 


0.808 


0.827 


30.7 


0.819 


0.821 


0.842 


44.7 


0.839 


0.841 


0.864 


52.8 


0.841 


0.843 


0.866 


62.5 


0.846 


0.848 


0.871 



3.5 Solution of Interior and Boundary Classes 

Consider the set of classes / U B = I w U I x U B w U B x . The boundary unitarity classes 
are 

B x = {j|a£ - an 2 - a^ 2 = 0, M x + 1 < J < M 2 } (3.65) 

and 

= { J\a J = 0, Mi + 1 < J < M 2 } (3.66) 

where M : = Floor[^/4/C(l + n/r 2 ) ], M 2 = Floor [y^] and C = -t/k 2 . The contri- 
bution to | Imr5| from the boundary unitarity class B can range from 0% to 100% and 
the contribution to |Imr 5 |, from the boundary unitarity class, can be selected without 
violating any of the constraints and this contribution can be made arbitrarily small. 

Consider the case with = 0.1E d , = 0.9E eb at ^/s = 52.8 GeV and t = 
—0.001 (GeV/c) 2 . The maximum contribution to | Imr 5 | is 34.7 where | Imr|| < 0.5 and 
| Imrf| < 34.2. The case with E^ = 0.01E cl , E^ = 0.99E cl , leads to | Imr 5 | < 11.6 where 
| Imrf| < 0.8 and | Imr^| < 10.6. The bound on | Imr^| falls when the fraction of E c i in 
the boundary unitarity class is reduced. The partial wave amplitudes in this region also 
become smaller in amplitude and contribute less to the bound on | Im r§ \ . The fraction 
of E e i in the boundary unitarity class can be further reduced until the contribution from 
this class to | Imrs| is negligible in comparison with the contribution from the interior 
unitarity class. In this limit the bound is, as before, 

|Imr 5 | <m^g~ f 36 ^^ _ iV x h (t) (3.67) 

\ ^tot / 
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or, under the approximation g ps of ot /(327rcr e i), 



where 



Imr 5 | < mwf x h(t) (3.68) 



(l + ^(2 + 9/8^)) 1/2 
/»(*) = . (3.69) 



The bound is identical to the bound when only the interior unitarity class is considered. 
A finite number of partial waves at low momentum transfer ensures a finite value for the 
sum 

lm4> 5 = J2J 2 a J 21 ( X -^ j2 ) ( 3 ' 7 °) 
and consequently a finite upper bound on | Im r 5 1 of less than unity. 

4 Error on Bound 

The upper bound on the imaginary single helicity-flip amplitude, modified by a kinemat- 
ical factor at zero momentum transfer, is given by 

i^WIFSH" (41) 



There are errors on all the experimental quantities in Eqn. (4.1) and consequently the 



upper bound on | Imr 5 | has an uncertainty. The experimental quantities g, cr to t and a cl 
have a nominal value plus an uncertainty: g±Ag, <7tot±Acrtot and <7 e i±A<7 e i. What are the 
values of Ag, Aa tot and A<r e i? Consider the value of Imr 5 at y/s = 52.8 GeV; | Imr5| < 
0.891 with g = 6.435 ± 0.14GeV~ 2 [I7|, a m = 42.906 mb JTJ and a cl = 7.407 mb JT| . 



The value of Ag is known but we must calculate Acr tot and A<r e i. 

A parameterization for the total and elastic cross section in elastic pp collisions [ |T8] ] 
allows a value for the cross sections to be found and a value of their uncertainties to be 
calculated. Each cross section is parameterized as 

a(p) =A + Bp n + C \og 2 (p) + D log(p) (4.2) 

where a is in mb and p is the laboratory momentum in GeV/c. The uncertainty in a is 
given by 
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1/2 



r AC) + \0D <AD) < (4 ' 3) 



The fitted parameters A, B, n, C and D are given in Table |j. 

Table 4: Fitted parameters for pp scattering 



Reaction 


A 


B 


n C 


D 


Ctot 


48.0 ±0.1 




0.522 ±0.005 


-4.51 ±0.05 




11.9 ±0.8 


26.9 ± 1.7 


-1.21 ±0.11 0.169 ±0.021 


-1.85 ±0.26 



Using Eqn. ( |4.3|) , and the values of the fitted parameters, the uncertainties Acr to t and 
Add can be written as 



A<7 to t = yO.01 + 2.5 x 10~ 5 log 4 (p) + 2.5 x 10" 3 log 2 (p) (4.4) 

and 

Aa el = ^0.64 + 2.89p- 2 - 42 ± 12.819p- 4 - 42 ± 4.41 x 10" 4 log 4 (p) ± 0.0676 log 2 (p) . (4.5) 

A laboratory beam momentum of p = k^fsjm = 1485 GeV/c at \/s = 52.8 GeV gives 
A(x tot = 0.463 mb, or 1.08% of a tot and Aa el = 2.345 mb, or 31.66% of a ch 

The uncertainty in Imrs is 
A Im r§ 



9 Im ^ 2 (A 9 y + ( ?^H) 2 ( Aa tot y + ( 9 -^) 2 ( Aa cl y . (4.6) 



At = 52.8 GeV, g = 6.435 ± 0.14GeV~ 2 , a tot = 42.906 ± 0.463 mb and a el = 7.407 ± 
2.345 mb. The uncertainty Almrs = 0.049 and the upper bound on | Imrs| is 0.891 ± 
0.049. 



The approximation 



can be used to write the bound on | Im r 5 1 as 



Imr 5 | < rnJ— (4i 
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at zero momentum transfer. The uncertainty A Im rs is simply 



/ (9Imr5 



) 



2 



A Im rs = 



(Ag) 2 . 



(4.9) 



At y/s = 52.8 GeV, g = 6.435 ± 0.14GeV" 2 , |Imr 5 | < 0.846 and AImr 5 = 0.009 or 
| Imr 5 | < 0.846 ± 0.009. The error on a c \ being large has relatively little effect on the 
uncertainty of the bound. 

5 Conclusion 

A bound on the size of the imaginary part of the single helicity-flip amplitude , in the CNI 
region, was obtained using the Lagrange variational method of Einhorn and Blankenbecler 
with o"tot, o"ci and diffraction slope expressed as equality constraints and unitarity expressed 
as two inequality constraints. An important feature associated with the bound is the 
number of partial waves at low momentum transfer. In the CNI region the number of 
waves is finite, that is, there is no singularity in \f^i as t — > 0. This feature ensures 
that a useful bound exists and, in fact, the bound limits | Imr 5 | to values less than unity 
near t = 0. With more constraints in the system an improved bound could be obtained 
although any additional constraints must be sufficiently well known experimentally and 
be computationally tractable. As the bound of 0.84 on the helicity flip amplitude ratio is 
less than (/i p — l)/2 = 0.896 at the high energies considered, the coefficient of (t c /t) in the 
expression for the asymmetry is constrained to be positive and therefore the analyzing 
power is positive for at least a significant part of the interference region. Though the bound 
is not sufficient to limit the polarization error to the recommended 5%, it does permit the 
use of proton proton elastic collisions in the CNI region as a relative polarimeter. 
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